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Abstract 

In the supersymmetric grand unified theories (SUSY GUTs), gauge bosons asso- 
ciated with the unified gauge group induce proton decay. We investigate the proton 
decay rate via the gauge bosons in the SUSY GUTs under the two situations; one is 
with extra vector-Hke multiplets, and the other is with heavy sfermions. It is found 
that the proton Hfetime is significantly reduced in the former case, while in the 
latter case it is slightly prolonged. Determination of the particle contents and their 
mass spectrum below the GUT scale is important to predict the proton lifetime. 
The proton decay searches have started to access to the 10^6 GeV scale. 



1 Introduction 



The grand unified theories (GUTs), which embed the Standard Model (SM) gauge group 
into a large single gauge group, are quite attractive, and so a variety of models of the 
theories are proposed since the earliest work based on the SU(5) symmetry group was 
presented by Georgi and Glashow in 1974 [1]. Among them, the supersymmetric grand 
unified theories (SUSY GUTs) are considered to be promising candidates since they realize 
the gauge coupling unification with great accuracy [2] as well as solving the hierarchy 
problem in the GUTs [S]!^. The supersymmetric version of the Georgi- Glashow SU(5) 
GUT, which is the simplest among the SUSY GUTs, is called the Minimal SUSY SU(5) 
GUT g]. 

One of the most distinctive features which GUTs predict is the existence of baryon- 
number violating interactions, such as proton decay. In the SUSY GUTs with R parity, 
for example, proton decay is induced by two different processes: the colored-Higgs and the 
X-boson exchanging processes. The colored Higgs triplets are introduced for the Higgs 
doublets in the Minimal Supersymmetric Standard Model (MSSM) to be incorporated 
into the SU(5) multiplets. The colored-Higgs exchange yields baryon-number violating 
dimension-five operators, which give rise to the dominant contribution to proton decay 
in the Minimal SUSY SU(5) GUTs [5]. On the other hand, X bosons are the gauge 
bosons in the unified SU(5) gauge group, and they induce proton decay through the 
dimension-six operators. Since there has been no experimental signal of proton decay so 
far, strong limits are imposed on these interactions, especially on the former one. Super- 
Kamiokande gives the bounds on the proton lifetime: r(p — > e"'"7r°) > 1.29 x 10^'^ yrs with 
219.7 kt-yr of data [6] and r(p K+iy) > 3.3 x lO^^ ^j^j^ ^72.8 kt-yr of data at 90% 
confidence level [7j. These bounds are so stringent that the Minimal SU(5) SUSY GUT is 
excluded since the colored-Higgs-boson exchange process yields too short lifetime in the 
p — )■ K^i? channel [SJ. This prediction is, however, quite mo del- dependent. In fact, many 
models are proposed in order to allow them to evade the experimental constraints. Such 
attempts [9] are accomplished by pushing the mass of the colored Higgs much heavier than 
the GUT scale and/or suppressing the dimension-five operators by a certain symmetry, 
e.g., the Peccei-Quinn (PQ) symmetry [10] in the four- dimensional GUTs. In the extra- 
dimensional models, the R symmetry is introduced to suppress the dimension-five proton 
decay [H]. With the extensions, the X-boson exchanging interactions become dominant 
and the proton lifetime turns out to be long enough to avoid the current limit. In this 
work, we assume such mechanism works and concentrate on the proton decay via the 
X-boson exchange. 

Nowadays, the weak-scale supersymmetry itself is severely constrained by the exper- 
iments at the Large Hadron Collider (LHC). Since there has been no signal of SUSY 
particles, the ATLAS and the CMS collaborations give severe limits on their masses, 
especially those of colored particles [12]. Moreover, the collaborations have recently an- 
nounced that they detected signatures of Higgs boson for the first time [13]. According 
to their results, the mass of Higgs boson is about 125 GeV. The result is welcome to the 
MSSM, which predicts the light Higgs boson, while the substantial radiative corrections 
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are required in order to raise the mass of Higgs boson up to 125 GeV. Since the corrections 
are increased as the SUSY particles become heavy, this situation imphes that the SUSY 
scale might be higher than the electroweak scale [H], and challenges natural solutions to 
the hierarchy problem. These results give constraints on models of the SUSY-breaking 
mediation. The low-energy gauge mediation model may be required to introduce several 
messenger multiplets so that the SUSY particles are heavy enough [H]. On the other 
hand, there have been several alternatives to explain the 125 GeV Higgs boson in an 
extension of the MSSM. Introducing vector-like supermultiplets to the MSSM is one of 
the simplest way to accomplish the purpose [16]. In this case, the quantum effects by 
the extra multiplets help to increase the Higgs-boson mass. With the multiplets being a 
representation of the grand unified group, the perturbative gauge coupling unification is 
still preserved with great accuracy. 

The MSSM with vector-like multiplets and high-scale SUSY models modifies the pro- 
ton decay rates. In this work, we examine it in detail. We will find that vector-like 
matters enhance the proton decay rate, while heavy sfermions extend its lifetime. This 
investigation indicates that the proton-decay experiments might provide a hint for the 
low-energy structure of the SUSY GUTs. 

We assume that the SUSY GUTs are realized in the four-dimensional spacetime, and 
not consider the extra-dimensional SUSY GUTs in which the gauge symmetry is broken 
by the boundary conditions. In the latter models, the main decay modes of proton depend 
on configuration of quarks and leptons in the extra-dimensional space [17j. However, the 
decrease / enhancement of the proton lifetime we show are almost universal even in those 
models. 

This paper is organized as follows. We begin by reviewing the SUSY SU(5) GUTs in 
Sec. m and discuss the proton decay via the X-boson exchange in the subsequent section. 
Then, we evaluate the lifetime of proton decay in the case of the MSSM with vector-like 
matters and heavy sfermions in Sec. S] and Sec. respectively. Section is devoted to 
conclusion. In Appendix, we present formulae for evaluating the long-distance part of the 
renormalization factors. 



2 SUSY SU(5) GUTs 



First, we review the SUSY SU(5) GUTs in order to illustrate the notation and conventions 
which we use in this article. Just like the Georgi-Glashow SU(5) model [T], we assume 
that the SM fermions, as well as their superpartners, are embedded in a 5 © 10 represen- 
tation. The multiplets $ and \E', which are the matter fields of 5 and 10 representations, 
respectively, are identified as 
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where all of the component fields are expressed in terms of the left-handed chiral super- 
fields, and primes show that the above fields are the gauge eigenstates. The super/sub 
scripts of U'^, and D'^ denote a color index with a = 1,2, 3. The chiral superfields 

E' and N' in $ as well as C/'° and D'"' in ^ form the SU(2)i doublets, respectively. 
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while D'^, and E' are the SU(2)l singlets. 

The SU(5) gauge theory contains the 24 gauge bosons and each of them corresponds 
to a component of a vector superfield, V^, with A = 1, . . . , 24 indicating the gauge index. 
By exploiting the fundamental representation of the SU(5) generators, T"^, we define a 
5x5 matrix of the vector superfields: V = V^T^. The components of the matrix are 
written as 
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where each component is expressed by the same symbol as that used for the corresponding 
gauge field. We collectively refer to X^, Y^, and their Hermitian conjugates as X and Y 
bosons, and use the following notation for them: 
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Here r, s, . . . denote the isospin indices. 

The Higgs superfields in the MSSM, on the other hand, are incorporated into a pair 
of fundamental and anti-fundamental fields as 
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(5) 



where the last two components are corresponding to the MSSM Higgs superfields. 
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and the new Higgs superfields Hq and Hca are called the Higgs color triplet superfields. 
The superpotential for the Yukawa couplings of quarks and leptons is given as 



W = -h'^eabcde 



(7) 
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where i,j = 1,2,3 indicate the generations and a,b,c,... represent the SU(5) indices. 
The Yukawa couphngs h^^ and /*-^ in Eq. ([7]) have redundant components and much of 
the degree of freedom is ehminated through the field re-definition of \1/ and $ [18]. We 
parametrize the couphngs according to Ref. [19] as 

= , (8) 

r = V:^fd, , (9) 

with Vij the Kobayashi-Maskawa matrix. The phase factors (pi are subject to a condition: 

V^i + + V^a = . (10) 

With the parameters, we express the matter fields in terms of the mass eigenstates as 
follows: 

U: = e-"^'U„ D[ = D,, E[ = V,,E,. (12) 

In the following discussion, we express interactions in the basis of mass eigenstates unless 
otherwise noted. 

Here, we do not take specific assumptions for symmetry breaking of SU(5) and the 
mass generation of the colored Higgs. They are related to suppression of the proton decay 
induced by dimension-five operators. 



3 Proton decay via the X and Y boson exchange 

Next we discuss the proton decay rate induced by the X and Y boson exchange. The 
couplings between the gauge bosons and the SM fermions are given as 

+ e-^'^^e"'^^(Q;)„,^^^Pi(nf)^ + h.c] , (13) 

where is the SU(5) gauge coupling constant and C indicates the charge conjugation. 
ers and e'^^"' denote the second and third rank totally antisymmetric tensors, respectively. 
This interaction Lagrangian causes proton decay. In the present case, the dominant decay 
mode is p — 7- Tr^e"*", and this decay process is induced by the following effective Lagrangian: 



Aff = -^e^'^'^aM [A^nh'^TPRd^e+PLu' + + \V^,\''){WrPLdPe+PRU^] . (14) 



where Mx is the mass of X and Y bosons and -Pl/r = ^(1 T 75)- The renormalization 
effects resulting from the anomalous dimensions of the operators are represented by 
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(2'] 

and , which we will evaluate below. The hadron matrix elements of the operators are 
evaluated by using the ordinary chiral Lagrangian methocQ [22j : 

{A^^MWTPRd^]PLU^\p{P. s)) = ^(1 + D + F)Plu,{p, s) , (15) 

and 

{Ae^ls4(^rPLd^]PRuMP, s)) = ^(1 + D + F)PnuM , (16) 

v2/,r 

where /^r is the pion decay constant, and an is defined by the equation 

(0|e„^^KCP^rf'^]Pz.M^b(p, s)) = anPLUpip, s) , (17) 

with C the charge conjugation matrix. The value of an is computed in Ref. [23] as 

aH(2 GeV) = -0.0112 ± 0.0012(3tat) ± 0.0022(syst) GeV^ (18) 

at the renormalization scale fi = 2 GeV. By using the matrix elements, we obtain the 
partial decay width r(j9 — )■ vr^e"'') induced by the effective Lagrangian in Eq. ( !T4|) : 

(19) 

where and m^r are the masses of proton and the neutral pion, respectively, and = 
g'i/4:7i . In the following calculation, we take = 130 MeV, D = 0.80, and F = 0.47 as 
in Ref. 

Now, in order to evaluate the proton decay rate, all we have to do is to determine 
the unified gauge coupling constant as, the X boson mass Mx, and the renormalization 
factors A^'* and A^^\ They are dependent on each GUT model. 

The unified gauge coupling constant is computed by solving the renormalization 
group equations (RGEs) for the gauge coupling constants ga (a = 1,2,3) in the SM gauge 
interactions. In this article, we consider the gauge coupling running up to the two-loop 
level, and exploit the DR renormalization scheme in order to respect the supersymmetry 
[21]. Furthermore, we adopt the definition of as as = g^^MQu-r) / 4:H with Mgut = 
1.5 X 10^6 GeV. 

In the MSSM, the two-loop renormalization group equations (RGEs) for the gauge 
coupling constants are given as 



dga^J_^,), gl 
^ dfi 167r2 " ^ (167r2)2 



i=t,b,T 



(20) 



^ Calculation of the matrix elements is also conducted by using the direct method [20] • Recent 
progress, in which the quenched approximation is not used, is reported in Ref. |21j . 
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where 
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and 



'26/5 14/5 18/5^ 
6 6 2 
4 4 



t, &, t) 



(22) 



with yi the Yukawa couphngs. Here, we use the SU(5) normahzation for the U(l) hy- 
percharge. Since the Yukawa couphngs enter into the two-loop level contributions to the 
gauge coupling RGEs, it is sufficient to consider the RGEs for the Yukawa couplings at 
one- loop level. They are given as 
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In the SM, on the other hand, the coefficients for the gauge coupling beta functions 
are [26] 
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The running of the Yukawa couplings in this case is given as follows: 
d 1 
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d 
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(26) 



Modifications in the coefficients in particular models are mentioned to in the following 
sections. 
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Now we deal with the renormahzation factors, A^^'^ and . The factors are expressed 
as the product of the long- and short-distance factors, i.e., 
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where Al and A^-* represent the long- and short-distance factors, respectively. The long- 
distance contribution Al is common to A^^^ and A^^\ and independent of the high-energy 
physics. Its value is evaluated as 

Al = 1.25 , (28) 

at two-loop level. Details of the calculation are given in Appendix. 

The short-distance factors are, on the other hand, model-dependent quantities. For 
instance, if there is no threshold between the electroweak and the GUT scales they are 
evaluated at one-loop level in Refs. [271 121] as follows: 
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while in the SM, 



73=3, 



73 



72 



(1)^ 11 
30' 



(2) 23 
=30 



"Try 



(2) 



72 



9 

4' 



7i 



(1) 



11 
20' 



7i 



(2) 



23 
20 



(29) 

(30) 
(31) 



The extension of the result to each case is referred to in the subsequent sections. 

For reference, we evaluate the proton-decay lifetime assuming the SUSY scale to be 
1 TeV, i.e., all the superparticles are assumed to have masses of 0{1) TeV. The result is 



t{p e+7r°) = 1.16 x 10 
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years . 



(32) 



Here we neglect the possible effects of the threshold corrections from particles whose 
masses are around the GUT scale. We also neglect them in the following calculations, 
since the effects are completely mo del- dependent. 



4 Proton decay with vector-like matters 

In this section, we discuss the grand unified models in which extra vector-like matters 
are added into the SUSY SU(5) GUTs. We assume that there exist 77,5 and nio pairs of 
chiral supermultiplets which transform as 5 -|- 5 and 10 -|- 10 representations, respectively. 
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and evaluate the proton decay rate for the cases. For brevity, all of the multiplets are 
assumed to have the same mass, M. As mentioned to in the Introduction, such models are 
motivated by the symptoms of Higgs boson with its relatively heavy mass, mh ~ 125 GeV, 
reported recently at the LHC [13]; if the additional multiplets with mass around weak 
scale couple to the MSSM Higgs bosons, the Higgs mass m/j can be raised on account of 
the quantum corrections [TB]. The gauge mediation requires the vector-like matters as 
the SUSY-breaking messenger. Thus, we take a wide parameter range for M. 

The existence of the extra vector-like multiplets modifies the RGEs for the gauge cou- 
pling constants, as well as the short- distance renormalization factors. The beta functions 
of the gauge coupling constants in the MSSM receive additional contributions [29]: 
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Here 6b^a^ and 6bf^ are the corrections to the leading and the next-to-leading order con- 
tributions for the beta functions of the gauge couphng constants, fei^^ and b^^fj , which are 
defined in Eq. fl?Il) . Here, we ignore the contribution from the Yukawa couplings of the 
vector-like matters at two-loop level for simplicity. From the one-loop contribution in 
Eq. (1331) . it is found that addition of vector-like matters equally changes the running of 
the gauge couplings, and thus maintains the gauge coupling unification. 

The corrections to the one- loop contribution, 6b^a \ also modify A^g^ and A^p . To be 
concrete, the renormalization effects of the energy scale above M are given as 
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These modifications result in alternations of the proton-decay lifetime. In order to 
parametrize the alternations, we define the following ratio: 
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r{p e+7rO)|„/o 



(35) 



where r(p — )• e~^7r°)|w/(w/o) represents the proton-decay lifetime with (without) vector-like 
matters. As seen form Eq. (fT9|) . this ratio does not depend on the X-boson mass but only 
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Figure 1: Ratio R as functions of vector-like matter mass. Left and right graphs show 
= 1, . . . , 5 and nio = 1, ... 3 in sohd hues from top to bottom, respectively. Light 
(dark) shaded region is excluded by the current experimental limit, r(p — > e'^vr'') > 1.29 x 
10^^ years at 90% confidence level |6], in the case of Mx = 1.0 x 10^^ GeV (2.0 x lO^*^ GeV). 
SUSY scale is set to be 1 TeV. 



on as, A^^^ and ^5^: 



«i[(4))V(AgV(i + |K.IT]„/ ■ ^''^ 

Here, w/ (w/o) again implies that the factor is for the case with (without) vector-like 
matters. 

In Fig. [H we plot the ratio R as functions of the masses of the vector-like multiplets. In 
the analysis, the SUSY scale is set to be 1 TeV. Each solid line in the left graph corresponds 
to the number of 5 -|- 5 multiplets = 1, 2, . . . , 5 from top to bottom without 10 + 10 
multiplets. In the right graph, on the other hand, the solid lines represent the cases with 
10 + 10 multiplets (nio = 1, 2, 3, from top to bottom) and = 0. The light (dark) shaded 
region is excluded by the current experimental limit, r(p — t- e+7r°) > 1.29 x 10^^ years 
at 90% confidence level [6] in the case of Mx = 1.0 x lO^^ GeV (2.0 x 10^^ GeV). On 
the whole parameter region in this figure, the unified gauge coupling constant, as, is 
less than 0.6, thus, perturbativity is maintained. In order to clarify the effects of the 
renormalization factors on the enhancement of proton decay rate, we also plot the ratio of 
the factors, [{A^^^f + (A^ni + \VM\/ /[(A^i^)' + (A^R^fil + \VM\/o, as functions 
of the masses of the vector-like multiplets in Fig. |2j Here we set ^s = 1, . . . , 5 from bottom 
to top and the SUSY scale to be 1 TeV. The renormalization factor in the absence of the 
vector-like multiplets is calculated as [(^1^'')^ + (y4^^)^(l + 1 V^rfp)^]„/o = 40 . These results 
indicate that the MSSM with vector-like matters might cause proton decay fast enough 



9 



o 

CO 
DC 



10 - 




lO'* 10^ 10^ 10^° 

Mass (GeV) 



10 



,12 



10 



,14 



Figure 2: Enhancement of renormalization factors as functions of vector-like matter mass. 
We take ns = 1, . . . , 5 from bottom to top. SUSY scale is set to be 1 TeV. 



to be detected by the current or future experiments. In particular, if there are three 5 + 5 
vector-like matters or is a 10 + 10 multiplet at the TeV scale, the SUSY SU(5) GUTs are 
excluded even by the present limit for Mx = 1.0 x 10^^ GeV. 

In addition, we calculate the ratio in the case where Ug pairs of 5+5+10-1-10 multiplets 
are introduced. This set of representations is the same as those of one generation of chiral 
matter fields and their corresponding complex representations. In this case, the beta 
functions are equal to those of ris = nio = Ug in Eq. f l33|) . After carrying out a similar 
calculation, we plot the ratio R against the masses of the multiplets in Fig. [31 Here, solid 
lines represent rig = 1, 2 multiplets of 5 + 5 + 10 + 10 from top to bottom. We again set 
the SUSY scale to be 1 TeV. From this figure, it is found that a pair of extra generations 
which have masses below 10^ GeV yield too large proton decay rate to be excluded by the 
current experimental limit in the case of Mx = 1.0 x 10^^ GeV. Although the threshold 
corrections around the GUT scale might alter the predicted values of proton lifetime, a 
growing tendency in the proton decay rate is independent of particular models. Therefore, 
the proton-decay experiments are extremely promising for constraining such scenario. 

5 Proton decay with heavy SUSY particles 

In turn, we discuss the SUSY GUT with sfermions having masses much larger than the 
TeV scale. Such a mass spectrum might be realized when the supersymmetry is broken via 
the anomaly- mediation mechanism [30], and there have been a lot of works which examine 
the scenario [31]. In this scenario, although a high degree of fine-tuning is inevitable, the 
SUSY flavor and CP problems are relaxed owing to heavy masses of sfermions |32] . The 
thermal relic scenario of dark matter in the Universe is still achieved [33], and the dark 
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Figure 3: Ratio R as functions of the masses of the multiplets. Sohd hues represent Ug = 
1, 2 multiplets of 5 + 5 + 10 + 10 representations from top to bottom. Light (dark) shaded 
region is excluded by the current experimental limit, r(p — )■ e"^7r°) > 1.29 x 10^^ years 
at 90% confidence level [6] in the case of Mx = 1.0 x lO^^ GeV (2.0 x lO^^ GeV). SUSY 
scale is set to be 1 TeV. 



matter is to be directly detected in the future experiments [3l]. In addition, this heavy 
scale SUSY scenario is also suggested by the recent LHC results since there has been no 
signal of superparticles and the ATLAS and CMS collaborations provide stringent limits 
on the masses of colored particles [12]. Moreover, in this case the SM-like Higgs boson 
mass easily reaches ~ 125 GeV as the radiative corrections are enhanced due to the heavy 
colored particles. 

In the following discussion, we assume that the masses of squarks and sleptons are 
around the scale of Msf, which is much higher then the electroweak scale, while the masses 
of gauginos, higgsinos, and Higgs bosons except the lightest one are C(l) TeV, which is 
denoted by Mgh hereafter. Therefore, between Mqh and Msf, the latter fields are to 
be added to the RGE analysis. Furthermore, we assume that their contributions to the 
gauge-coupling running through the Yukawa couplings, as well as those with sfermions 
running in the loops, are negligible. Hence, the beta functions of the gauge couplings in 
this region are given as 

/23/5\ /104/25 18/5 44/5\ 

6^ = -1 , 6l? = 6/5 32 12 , (37) 
\ -5 / \ 11/10 9/2 22 / 

with Cat the same as the SM ones. Since each one-loop contribution b^a'^ in the above 
equation differ from that in Eq. (!2T|) by the same number (2 in this case), the perturbative 
gauge coupling unification is again preserved in the present case. 
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Figure 4: Proton lifetime against Msp- Solid, dashed, and dotted lines correspond to the 
cases where Mgh is set to be 1, 3, and 10 TeV, respectively. X-boson mass is taken to be 
1.0 X 10^^ GeV. 



In Fig. m we plot the lifetime of proton against Msf- Solid, dashed, and dotted lines 
correspond to the cases where Mgh is set to be 1, 3, and 10 TeV, respectively. The 
X-boson mass is taken to be 1.0 x 10^^ GeV in this figure. It is found that the proton- 
decay lifetime is slightly extended , although the enhancement factor is less than two. 
The contribution of the renormalization factors to the alternation of the proton lifetime 
is less significant in this case than that in the case discussed in Sec. HI It is at most a 
few % on the whole parameter region in this figure. Compared with the case discussed 
in the previous section, the present situation does not so much change the proton decay 
rate. Hence, searching for the proton decay is still stimulating even for the heavy SUSY 
scenario. 

6 Conclusion 

We have studied the proton decay rate under the two situations; one is the MSSM with 
vector-like multiplets, and the other is the MSSM with heavy sfermions. It is found that 
the proton lifetime is significantly reduced in the former case, while in the latter case 
it is slightly prolonged. In any case, the proton-decay experiments are, together with 
the LHC experiment and other precision measurements, expected to shed light on the 
supersymmetric grand unified models, as well as the supersymmetry itself. 
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Appendix: Long-distance part of the renormalization 
factors 

In this Appendix, we demonstrate the evaluation of the long-distance contribution to 
the renormalization factors, e.g., Al in Eq. (127|) . First, we write down the two-loop 
renormalization group equations for the strong coupling constant below the electroweak 
scale: 



dgs 



with Qs the strong coupling constant and 

2 



(167r2)' 



;b29. 



5 

s 1 



hi 



11 



-N. 



102 



38 



f ' 



(38) 



(39) 



where Nf denotes the number of quark flavors in an effective theory. 

The long-distance factor, A^, is determined by the ratio of the coefficients for the 
effective operators at the scale of mz and 2 GeV: 

C{2 GeV) 



with the coefficient C(/i) satisfying the following RGE at two- loop level [35] 



4c(.) = 

The solution of the equation is 



as /4 4, \ a, 
4— + - + -Nf — — 



c(/i) 




bl 


C(/io) 







Airbi + b2as{jj,) 
A-nhi + 62as(/io) 
with hi and 62 given in Eq. (!39|) . Thus, Ai'is given as 

a,(2 GeV) + 
a,(mfe) + ^ 



6+2JVj 
962 



Ar 



'as{2 GeV)" 


6 

25 




6 

23 


as{mb) 




_as{mz)_ 





«s(^b) + ^ 

as{mz) + ^ 



(40) 



(41) 



(42) 



(43) 



and numerically it turns out to be 



Aj 



1.25 . 



(44) 
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